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Abstract. In this paper, we present the implementation of Bloch oscillations in an atomic interferometer 
to increase the separation of the two interfering paths. A numerical model, in very good agreement with 
the experiment, is developed. The contrast of the interferometer and its sensitivity to phase fluctuations 
and to intensity fluctuations are also calculated. We demonstrate that the sensitivity to phase fluctuations 
can be significantly reduced by using a suitable arrangement of Bloch oscillations pulses. 
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1 Introduction 

Atom interferometry is nowadays a corner stone in high 
precision measurement. It has been used to determine the 
fine structure constant |1I2) , as well as the Newton's grav- 
itational constant |3I4) , to test general relativity [5] and to 
measure gravity [617] . In these applications, the interfer- 
ometers rely on momentum exchange between light and 
atoms to split and recombine atomic wave-packets. The 
sensitivity of the measurements is proportional to the spa- 
tial separation of the wave-packets in the interferometer. 
As the beamsplitters arc performed with two photon Ra- 
man transition, the velocity difference between the arms 
of the interferometer Sv is equal to 2v r (where v r = hk/m 
is the recoil velocity, k is the wavenumber of light and m 
the mass of the atoms). The sensitivity of the interferom- 
eter is thus proportional to 2v r T where T is the typical 
duration of the interferometeiQ 

Until recently, the main way to improve the sensitiv- 
ity of interferometers in order to achieve higher precision 
measurements was to increase the duration T . However, 
it requires a careful control of phase noise of the lasers as 
well as vibrations in order to keep a good signal to noise 
ratio. Transversal spread of the atomic cloud or motion 
of the cloud inside the science chamber set a limit to this 
interaction time. Interaction times of the order of several 
hundreds of ms are achieved [B] . 

Another method to increase the separation of the wave- 
packets and so the sensitivity of the interferometer con- 
sists in using a large momentum transfer pulse with a 
separation Sv > 2v r . High order Bragg diffraction of mat- 



For a gravimeter, the velocity change due to gravity is also 
proportional to T, which leads to a sensitivity to gravity scaling 
asT 2 
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ter wave can be used to increase Sv [5]. However, the laser 
power required for Bragg diffraction increases sharply with 
Sv |5]. Recently, a double-diffraction technique has been 
proposed to enhance the area of a Raman atom interfer- 
ometer [TU]. Increasing Sv can be also done with Bloch 
oscillations. It has been suggested and demonstrated by 
the group of W.D. Phillips on a Bose Einstein conden- 
sate [TT]. The coherence of the acceleration with Bloch 
oscillations (BO) is a key point of the process [IS]. More 
recently, a group at Stanford/Berkeley and our group have 
implemented this method by inserting the so called large 
momentum transfer beam splitter (LMTBS) in an atom 
interferometer [1319] . 

The large momentum transfer beam splitter (LMTBS) 
is realized with a suitable combination of two processes. 
The first one, realized with usual methods (Bragg pulse 
in the group of Stanford/Berkeley and Raman pulses in 
ours), implements a separation between two wavepackets. 
The second process is used to selectively and coherently 
accelerated only one of the wavepackets with Bloch oscil- 
lations in order to increase Sv. This separation keeps the 
coherence between the two wave packets. 

Our experiment where we demonstrate the realisation 
of a LMTBS has been described in a letter [13]. In this 
paper, we present a model that help us to calculate pre- 
cisely the efficiency of the LMTBS (first part) and also the 
phase of the interferometer (second part). The calculation 
of this phase is a key point to understand the sensitivity 
of the interferometer (third part). 



2 Bloch oscillations 

In this section we describe the Bloch oscillation process 
that is used to increase the separation of atoms in the 
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interferometer. The phenomenon of Bloch oscillation is a 
well understand tool to coherently accelerate ultra-cold 
atoms[T3]. In order to realise a LMTBS, the difficulty is 
not only to efficiently accelerate atoms but to be able to 
accelerate only atoms from on arm of the interferometer, 
keeping the coherence between the two arms. 



2.1 Phenomenological description of the acceleration 

The Bloch oscillation acceleration is based on coherent 
transfer of momentum between an accelerated optical lat- 
tice and atoms. The optical lattice results from interfer- 
ence of the two lasers beams of the same intensity /, de- 
tuned by A from the one photon transition. The amplitude 
of the resulting standing wave is 41. Atoms endure a pe- 
riodic potential U(x) due to light-shift (a.c. stark effect) 
which can be expressed as: 



' V # /V 



U(x) 



U 



(f + cos(2fca;)), 



(1) 



where k = 2tt/X = 2ttv/c is the wavevector of the laser, 
Uq = h-^j- with r representing the natural linewidth 
of the excited state and Is its saturation intensity. In the 
following, the constant part Uq/2 which does not play any 
role for the explanations is removed from every equation. 
Because of the periodicity of the potential, the eigenener- 
gies of the system show a band structure [T5] . Each eigen- 
statc is described by a band index n and a quasimomen- 
tum q defined modulo 2hk [IB]. The Fig. [T] represents this 
band structure, where the first Brillouin zone is unfolded. 
The splittings between bands at the center (q — 0) and at 
the edge (q — hk) of the first Brillouin zone occurs where 
atoms are resonantly coupled by U{x). 

In a quantum optics view, the avoided crossing (Bragg 
resonance) between the first and second bands is con- 
nected to the resonance of the two photon transition be- 
tween states of momentum p — ±hk. Similarly, the split- 
ting between the second and third bands is due to a reso- 
nant four photon transition between states of momentum 
p = ±2hk, and so on for the n th and (n + I ) th band. The 
order of the transition grows with the band index, thus 
the higher the band index is, the weaker is the splitting 
between the bands. 

When the atoms in the first band are subjected to 
a constant and uniform force, their quasimomentum in- 
creases linearly with time. If the force is weak enough to 
induce adiabatic transition, the atoms stay in the first 
band and therefore have a periodic motion (Bloch oscil- 
lations [H]). The time necessary for a variation of the 
quasimomentum of 2hk is the period of the oscillations. 
For atoms in higher band, for a constant acceleration, the 
probability to make an adiabatic transition (stay in the 
same band) will be higher for low values of the band in- 
dex. 

Bloch oscillations were initially introduced in the case 
of a lattice and a uniform force. In our experiment, a 
constant acceleration is applied to the lattice. This ac- 
celeration acts like a force in the frame of the lattice (see 
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Fig. 1. Band structure of the optical lattice. This graph repre- 
sents the band energy in the frame of the lattice as a function 
of the quasi-momentum. Because quasi-momentum is denned 
modulo 2hk, one can think of this graph as an unrolled cylinder 
where the vertical dot-dashed lines are overlapping. 



Ref. [H]), and an equivalent of Bloch oscillations can be 
observed [14] : the adiabatic transition corresponds to a 
change of velocity of 2 recoils, whereas a non-adiabatic 
transition does not change the velocity of the atom (free 
particle). 

In our LMTBS, one wavepacket is loaded in the first 
band (A, see Fig.rjJ and the second in the third band (B). 
The acceleration is small enough so that atoms in the first 
band are likely to make an adiabatic transition, but strong 
enough so that atoms initially in the third band change 
band. The atoms in the first band are periodically acceler- 
ated by 2 recoils momentum. In the same time, atoms in 
the third band are not accelerated. This process allows us 
to selectively accelerate atoms depending on their initial 
velocity. 

Note that in our experiment, atoms with different ve- 
locities have also different internal states because of the 
Raman transition used to create the initial splitting. The 
amplitude of the lattice is then different for the two wavepack- 
ets. This effect could be used to select the wavepacket 
which is going to be accelerated. However for our experi- 
mental parameters, the difference in light shift is negligible 
and so the selection of the accelerated wavepacket relies 
only on its initial velocity and not its internal state. 



2.2 Theoretical description of LMTBS 

In order to optimize the efficiency of the LMT beamsplit- 
ter, the evolution of the atomic system has been com- 
puted. In the laboratory frame, the Hamiltonian of the 
particle is: 

H=^- + V(x,t) (2) 



where V(x, t) is the accelerated periodic potential. If we 
denote by Uo(t) the amplitude of the lattice and by Sv(t) 
the frequency difference between the two laser beam. Then 
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v(t) — 2k5v(t) is the velocity of the lattice and 



V(x,t) = 



U (t) 



cos 2fc 



x- J v(t')dA\ (3) 



We should note that, even if the lattice is accelerated, 
in the laboratory frame, this Hamiltonian is still periodic 
in space. Therefore the quasimomentum is still a good 
quantum number and is conserved. This Hamiltonian can- 
not be solved efficiently, because of the fast time variation 
of the potential. One could write the Hamiltonian in the 
accelerated frame. This Hamiltonian will be time inde- 
pendent, but will no longer be periodic, and the quasi- 
momentum will not be conserved. 

We use another and convenient transformation which 
consists in translating the wave function in position space 
by the quantity X = J v(t')dt': 

x =x + X(t) 
The Hamiltonian becomes: 



H(t) = 



(p-mv(t)f Uo(t) 

2m 2 



cos (2kx) , 



(4) 
(5) 



(6) 



where the energy shift mv(t) 2 /2 is left out. This Hamil- 
tonian described the Bloch oscillation in the "solid-state 
physics" point of view [H] . The advantage of this presen- 
tation is that the Hamiltonian is time independent when 
the lattice is moving at a constant velocity. 

If we start with a plane wave of momentum po (and 
therefore of quasi- momentum qo = po modulo 2hk)), the 
wavefunction will remain in the subspace of states of quasi- 
momentum qo. We therefore choose to solve the problem 
in the basis of the plane wave of momentum po + 2lhk, 
I e Z which is a basis of this subspace. 

To solve numerically the problem i npractice, the mo- 
mentum basis is restricted to a finite number of I, \l\ < 
iV cu t. The value of iV cut is set by the maximal velocity 
of atoms in the optical lattice. It results from the sum 
of the number of oscillations (velocity of the atom) and 
the width of the velocity distribution of the Bloch states 
(Ap/hk). This width Ap is set by the Wannier function in 
momentum space [18 which is given by the wave function 
of an atom trapped in a single well without tunnelling be- 
tween sites. The frequency uj of the trap is proportional 
to y/\Ua\E r . For the ground state of an harmonic oscilla- 
tor, Ap = yJhmuj/2. We obtain that Ap/hk ~ {/U /E r 
|18ll9j . In our calculation, we have used iV cut = 10 for 2 
oscillations and Uq ~ 8E r . 

The Hamiltonian is calculated using the dimensionless 
parameters k = Uo/8E r (where E r — ^-^- is the recoil 
energy) and p(t) = (po + mv(t))/hk. All the units are 
scaled so that h = k = m = 1. 



The expression of the Hamiltonian H and the ket \tp) 
are reduced to : 
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C-l 



C() 



C-l 



(7) 



where p\ = p(t) + 21. Each coefficient q represents the 
amplitude of probability for the state to be in the plane 
wave of momentum 2lhk. 

The LMT pulse is switched on (and off) adiabatically 
to maximize the number of atoms transferred from an 
atomic state with well known momentum po (plane wave) 
to a Bloch state with a quasimomentum qo — pa in the 
band n (and reciprocally) [18 . In between, the lattice is 
accelerated to transfer a given number of recoils to one 
component of the wavefunction. We denote by t a diab the 
duration of the linear ramp used to switch on and off the 
lattice, T acc the duration of the acceleration, 2N the num- 
ber of transferred recoils, Uq the maximal amplitude of the 
lattice and po the initial momentum of the atoms. 



2.3 Optimization of the acceleration 

The optimization of the acceleration results from a com- 
promise between an adiabatic acceleration of the atoms in 
the first band and a strong enough (non-adiabatic) accel- 
eration of the atoms in the third band. Reciprocally, for 
a given acceleration, the amplitude of the lattice should 
be in a given range in order to satisfy this latest compro- 
mise. In our experiment, it is easier to vary the amplitude 
of the lattice. Therefore, we have plotted on the Fig. [5] 
the probability for an atom to stay in its band as a func- 
tion of the lattice amplitude n for the first (solid line) and 
third (dashed line) bands. There is clearly an intermediate 
regime where the probability rjn for an atom to stay in 
the first band is high whereas the probability 7734 for an 
atom to leave the third band (and reach the fourth one) 
is also high. The efficiency of the LMT is calculated as 



V = VuVte- 



(8) 



It sets the contrast of the interference pattern of the two 
wavepackets. The value of r\ is plotted on Fig. [5J For our 
experimental parameters (an acceleration of 2 recoils in 
200 /is), the maximum efficiency is about 98% for k — 
1. This maximum 77 maX j and the optimal value K opt of k 
depends upon T acc . The values of rj mayL and K op t versus 
T acc are plotted on Fig.[3l For an infinite duration of T acc , 
the efficiency is closed to one while the value for K opt is 
nearby zero. 
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Fig. 2. Transfer probability as a function of the maximal op- 
tical depth k, = Uo/8E r of the lattice. The acceleration is in 
200 /is for N=2 recoils. Solid line: transfer probability for the 
first band 7711 ; dashed line: for the third band, 7733 ss 1 — 7734; 
diamond: Efficiency of the LMT pulse, 77 — 77117734. 



This behavior is understandable using the Landau- 
Zener criteria which can be used to calculate transition 
probabilities in the weak binding limit (k -C 1), [T2l : 



?/Tl 

7?34 



1 — exp 



exp 



768a 



(9) 

(10) 



where a = m 2 /h 2 k 3 dv/dt, is the dimension-less acceler- 
ation of the lattice. We clearly see that for a small value 

6 2 

of k, there is a value of a so that =g=- <C 1 and — 3> 1 

] _ t 768a a 

leading to efficiencies 7711 and 7734 close to one. In fact, 
there is no theoretical limitation to the weakness of k. 
There is always a suitable value of a which sets 77 close 
to unity. However, this suitable value of a becomes very 
small, which means that an infinite amount of time is re- 
quired for the acceleration. In an other way, with a given 
value of a, the higher is a, the lower is 77 with an optimized 
/-r. 

In the experimental point of view, the duration of LMT 
pulse can not be extend as much as we want. This dura- 
tion has to be smaller than the time separation between 
the 7r/2 pulses of our interferometer. The time separa- 
tion between n/2 pulses results from a compromise be- 
tween the narrowest width of the atomic fringes (and so 
the resolution of the interferometer) and the vibrations 
of the experimental set up. Practically, the comparison of 
the theoretical and experimental results presented in this 
paper is done with T acc = 0.2 ms (for two oscillations). 
We will briefly explain at the end of this paper a way to 
transfer a larger number of recoils using a non constant 
acceleration. 



2.4 Adiabatic loading of atoms in the lattice 

In the previous section, the efficiency 77 of the separation 
has been calculated assuming that atoms are initially in a 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Tacc [ms] 

Fig. 3. Efficiency of the LMT pulse as a function of the total 
duration of the acceleration for two Bloch oscillations. For each 
value of the duration, we calculate the value K opt of the ampli- 
tude k of the lattice that optimizes the efficiency. This value is 
plotted in the inset of the figure. The "saw tooth" shape of K op t 
(and therefore of Tj op t ) is due to small oscillations between the 
bands. They do not appear in the Landau-Zener model which 
doesn't take into this effect (lost atoms cannot go back into 
their initial band). 



well defined Bloch state. Yet, in our experiment, we deal 
with a subrecoil velocity distribution, thus the best de- 
scription of the initial state is in a plane wave basis. To 
load the atoms into a Bloch state, the lattice depth is in- 
creased adiabatically. Atoms in an initial momentum po 
are loaded into the Bloch state \n,po), where the band 
index n is the nearest integer above \po/Hk\. 

The amplitude of probability rj t to load atoms in a 
given band is defined as \(ip\ n,Po)\ where \ij>) is the wave 
function and |n,po) the Bloch wave function. This proba- 
bility depends upon the loading time £ a diat>, the final depth 
of the lattice E/rj (the amplitude of the lattice is linearly 
increased) and the momentum po- The Fig. 2] shows, as 
functions of the duration £ a diab, the probabilities 774 to 
load atoms either in the first or in the third band. With 
an initial quasi- momentum go = 0.5, the probability rj t is 
higher for the third band than for the first one because the 
coupling with the lattice gets weaker with the increase of 
band index. 

On Fig.[5l this efficiency 774 is plotted versus the initial 
quasi-momentum for i a diab = 150 /as and for Uq — 8>E r . 

In the case of atoms loaded in the first band, the effi- 
ciency is maximum at the center of the first Brillouin zone 
(po = 0). At the edge of the first Brillouin zone (po = hk), 
rjt is exactly 50%. In this situation, it is indeed not possible 
to adiabatically transfer the atoms into the optical lattice, 
because the initial state is degenerate. The same phenom- 
ena occurs on the two edges of the third zone (po = 2hk 
and po = 3hk). The same difficulty happens at the end 
of the LMT pulse when the atoms are transferred back 
from the lattice to a well defined momentum state (plane 
wave). 
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Fig. 4. Transfer probability for the loading of atoms in the 
optical lattice as a function of loading time iadiab. The quasi- 
momentum is go = 0.5 and the final amplitude of the lattice is 

Uq = SE r . 
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Fig. 5. Transfer probability for the loading of atoms in the 
optical lattice as a function of the initial quasi-momentum go 
in unit of hk for atoms in the first and third band of the lattice, 

iadiab = 150 /IS and K = 1. 



2.5 Total efficiency of the LMT 

We have calculated the total efficiency r] to t with the fol- 
lowing parameters: i a <jiab = 150 fis, T acc = 200 /is, N = 2, 
Uq = 8E r . Those parameters correspond to an optimiza- 
tion of the efficiency, keeping the total time (defined as 
2i a diab + T acc ) constant but varying the amplitude of the 
lattice and i a diab (see Fig. E]). 

The total efficiency is plotted as a function of the ini- 
tial momentum po on Fig. [7] It is computed including the 
contributions of the adiabatic loading (and unloading) of 
the atoms in the lattice. In the first band, the state of 
initial quasimomentum qo = I and in the third bands, 
the states of initial quasi-momentum go = and go = 1 
are degenerated states. Therefore these atoms cannot be 
adiabatically loaded in the lattice. Consequently, it is im- 
portant to optimize the width of the initial momentum 
distribution such that it fits in the range where the process 
is very efficient. However, the wider is the initial velocity 
distribution loaded into the LMT pulse, the higher is the 
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Fig. 6. Time evolution of the amplitude k = Uo/8E r of the 
lattice and of the velocity v of the lattice for a beam splitter 
with 2 oscillations. These parameters correspond to an accel- 
eration of about 100 m.s -1 . 
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number of atoms that contributes to the interferometer 
and so the larger is the signal to noise ratio. Nevertheless, 
note that the total efficiency is larger than 95% on a wide 
range of po ■ 



3 Phase shift of the interferometer 

In the preceding section, we have calculated the efficiency 
of the LMT beamsplitter. This efficiency is one of the key 
parameter for calculating the contrast of the interferome- 
ter. Another parameter of the interferometer, is the phase 
shift. In the following section, we describe in detail how it 
can be calculated. 



3.1 Ramsey-Borde interferometer 

Let us briefly recall the principle of a Ramsey Borde in- 
terferometer realized with Raman pulses (see Fig. |5|). It 
is based on a succession of four 7r/2 Raman pulses (or, 
a succession tt/2-tt-tt/2). Each Raman pulse couples two 
states (labeled \gi) and \g2}) of the hyperfine structure of 
the ground state of the atom via a two photon transition 
in a A scheme. The two laser used for the transition are 
counterpropagating, therefore, the two states have differ- 
ent velocities. This allow us to split the initial wave packet 
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in two. The relative phase between the two states is deter- 
mined by the phase of the laser beams used for the Raman 
transition. The second and third pulse are used as mirror 
to modify the trajectory of atoms and the fourth pulse is 
used to recombine the two wavepackets. 

A simple way to precisely understand the output fringes 
pattern of the interferometer is to consider the contribu- 
tion of each pair of 7r/2 pulses. After the first pair of 7r/2- 
pulses, the velocity distribution of atoms in state \gi) fol- 
lows a Ramsey fringes pattern (see Fig. [5]). Thanks to 
the Doppler effect, the frequency difference between Ra- 
man beams sets the center of the Ramsey fringes pattern 
and select the most populated velocity class in the initial 
velocity distribution. Atoms which have not been trans- 
ferred are eliminated. The second set of 7r/2 pulses, also 
Doppler sensitive, is used to probe the selected velocity 
distribution. The resulting signal of the interferometer is 
the product of these two Ramsey distributions. By scan- 
ning the frequency difference of the lasers used for the 
second set of pulses, the second Ramsey pattern is trans- 
lated in velocity space. The resulting signal versus laser 
frequency difference is the convolution of the two distri- 
butions created with each pair of 7r/2 pulses. The center 
of this convolution profile gives access to the velocity of 
the selected atoms. 

The resolution of the velocity measurement is given 
by the fringe width. For a counterpropagating Raman 
transition, the Doppler shift of an atom at velocity v is 
2kv. The frequency splitting of the fringes in the convo- 
lution profile is 27r/T Ramsoy , i.e. A/2T Ramscy in velocity 
space. This expression can also be understood as resulting 
from the interference of two selected wave packets: they 
are spatially separated by AX — 2ivX Ram sey In velocity 
space, the interference produces fringes with a separation 



tAX 



2T R 



There are different ways of modelling an atomic inter- 
ferometer. One can, for example, use the Feynman path 
integral method [20121) . However, this method does not fit 
very well with the description of Bloch oscillation made in 
the first section of this paper, where we use plane wave. 

We therefore describe the atom interferometer using 
also a plane wave basis. We will first introduce how this 
description can be used for a regular Ramsey-Borde inter- 
ferometer and then extend it to the LMT interferometer. 



3.2 Modelisation of the interferometer 



We are going to calculate the evolution of a wave packet 
through the interferometer. This wavepacket will be de- 
scibed as a superposition of plane waves with a mean mo- 
mentum pq. During the interferometer the packet will be 
split into two paths (labeled A and B). Looking only to 
the evolution of a plane wave one can calculate the in- 
terference fringes and also the positions of the center of 
the wave-packet, given by the derivative of the phase with 
respect to the momentum po- 
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Fig. 8. Scheme of the usual Ramsey Borde interferometer. A 
first set of two tt/2 pulses, separated by a delay TR amsoy , se- 
lect a Ramsey fringe pattern in velocity space (1). A second 
set of two tt/2 pulses, with the same delay is used to close 
the interferometer. The probability for an atom to be trans- 
ferred by the second set of tt/2 pulses from one state to the 
other is recorded as a function of the difference between the 
frequency of the laser during selection 5 sc i and measurement 
<5 mC as (2). This interferometer measure then the velocity vari- 
ation between the selection and measurement pulses. 
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Fig. 9. Scheme of an interferometer using LMT beamsplitters. 
With this scheme four LMT Bloch pulses (labelled B) are in- 
serted between the Raman tt/2 pulses. In the experiment, we 
use a more sophisticated scheme described Fig. [10| 



The output of the interferometer can be precisely cal- 
culated. We have to solve the Schrodinger equation: 



in d -f = H { t m 



(ii) 



where the Hamiltonian H(t) is the sum of the free Hamil- 
tonian Hq and the Raman interaction Vr(£) with 



H = \-huj 12 \g 2 )(g2\ 

Am 



(12) 
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where huJi 2 is the energy difference between the two states 
\g 2 ) and \g x ) and 



Vk(t) = 



Q{t) 



j(4>B.(t)+K H x) 



(g\ + h. 



(13) 



where Q(t) is equal to the Rabi pulsation of the two pho- 
ton Raman transition (j?o during the pulses and else- 
where), hk e fi is the effective recoil transferred by the Ra- 
man pulses and <fin(t) is the phase difference between the 
two lasers. 

We choose to solve this Hamiltonian using the plane 
wave basis. We then have a two levels system \gi,po) and 
\g2,Po + Kk c s). We also write the phase of the laser as 
cf) R (t) = uj 12 t + 0(t) and change the ket \g 2 ) to e 4 " 12 ' \g 2 ) 
(rotating frame). The Hamiltonian is then this 2x2 matrix: 



H 



fl(t) e -**W/2 bo + hk cS ) 2 /2m 



(14) 



This Hamiltonian contains the phase shift induced by the 
kinetic energy and the phase transferred from light to 
atoms during a Raman transition. Assuming that <p(t) is 
almost constant during the pulse, we obtain that from \g) 
to |e) we add a phase 0(t) and from |e) to \g) we remove 
this phase. 

The final proportion of atoms transferred to state |e) 
will vary as sin(<£) with <P — 4>A~<t>B is the phase difference 
of the two paths with: 

0A = 01 - 03 + <f>f(Po) + ^ OS (P0) (15) 

<pB = 02 " 04 + S B O1 (PO) + 0r S (Po) (16) 



where <j>i — 0(£j) are the phases at the time £1,2,3,4 of the 
four Raman pulses and 0^/^ es (po) are the phases due to 
kinetic energy: 



0a'(po) 
0fbo) 



. Ill OS 

9 A 



(Po) 

bo) 



r. 



Ramsey ~ 



(p + hk e 



2mh 



= T, 



Ramsey 



2mh 



The phase can also be written as 



$ = ^lascr + $ 



kinetic 



(17) 
(18) 

(19) 



with0 lascr = (0i-03 )- (02 -04) and kinetic = (4>a(Po) + 

0r s (po))-(0 s B cl bo) + 0g cs bo)). 

We observe that ^kinetic = when the velocity of 
atoms is constant between the pulses. In the case where 
there are uniform forces that change po, we have to inte- 
grate the kinetic energy and we obtain that: 

^kinetic = — / \po{h + t) - p (t 3 + r))dr (20) 



The phase is a direct measurement of the momentum vari- 
ation between the first pair of pulses and the second one. 
Let us consider the case where there is a small change 
of velocity Sv between the second and third pulse. The sen- 
sitivity of the interferometer is proportional to the deriva- 
tive of <P with respect to 5v. 



Using the more general value of ^, and the fact that 
0a'(po) = 0r s (Po) and s B cl (p o ) = 4T s (Po), we calculate 
that 



dSv 



<90T 90f 



dp dp 



(21) 



This equation tell us that in order to compute the sen- 
sitivity of the interferometer, we need only to calculate 
precisely the phase shift accumulated during half of the 
interferometer. Provided that the forces are uniform, this 
phase can be calculated for plane waves only. 

Further more, in quantum mechanics, the mean posi- 
tion (X) of a wave-packet can be calculated as the deriva- 
tive of the phase with respect to momentum: 



<*>-*i 



(22) 



The equation [5T] shows directly that the sensitivity is pro- 
portional to the spatial separation of the wave-packet: 

— = 1 -((X A )-(X B )) (23) 



4 Large area interferometer 

The interferometer based on large momentum transfer 
(LMT) Bloch pulses is similar to the usual Ramsey Bordc 
interferometer. Raman pulses are used in the same con- 
figuration and the LMT Bloch pulses are simply inserted 
between the two 7r/2 pulses used either for the selection or 
the measurement of the velocity. A first implementation 
consists in accelerating and then decelerating the atoms of 
one arm of the interferometer for the selection (and then 
the atoms of the other arm for the measurement). The 
temporal sequence of such an interferometer is presented 
on Fig. M 

In the following section, we will present how this ac- 
celeration is implemented, and then the numerical calcu- 
lation of the phase shift <P of such an interferometer. 



4.1 Implementation of LMTBS in the set-up 

In our experiment, another arrangement slightly more com- 
plicated than the one on Fig. [HI is used: we accelerate suc- 
cessively each arm of the interferometer with two LMT 
pulses of opposite directions and then decelerate them 
with two other pulses. As we will see later, in such a way 
the contribution of some systematic effects is reduced be- 
cause the two arms are illuminated with Bloch beams in 
a more symmetric way. The whole interferometer used in 
our experiment is presented on Fig. 1101 It is realized with 
four Raman tt/2 pulses and eight LMT Bloch pulses (la- 
beled B on the figure [TU|) : each 7r/2 pulse is associated 
with two LMT pulses used to accelerate selectively each 
arm of the interferometer. 

After the first 7r/2 pulse, atoms arc in a superposition 
of states \gi,po) and \g 2 ,po + 2hk). A first LMT pulse is 
used to accelerate atoms with initial momentum pq + 2hk. 
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Fig. 10. Temporal sequence of the interferometer for N=2. 
From bottom to top: intensity of the Raman (dashed line) and 
Bloch (solid line) beams; velocity of the atoms in the first and 
second arms of the interferometer; trajectories of atoms in the 
two arms. 



Their momentum became po + 2hk + 2NHk after N os- 
cillations. Then atoms with initial momentum po are ac- 
celerated in the opposite direction. Therefore a superpo- 
sition of |g2,Po + 2S.fc + 2Nhk) and |<7i,po — 2Nhk) states 
is created. After a delay TR amsey , the same sequence of 
two LMT pulses is repeated but in reverse order, to bring 
back atom states in the initial superposition of \gi,po) 
and |<72> Pa + 2hk). Finally, the sequence is ended with a 
second Raman 7r/2-pulse. The wavefunction of atoms in 
the state I32) results from the interference between atoms 
transferred during the first 7r/2 pulse and atoms trans- 
ferred during the second one of the selection sequence. 

The interferometer is then closed in position space us- 
ing a second set of Raman and LMT pulses identical to the 
first one. As for the regular interferometer, the resulting 
signal of the interferometer (population in \g\) and I52)) 
vary as sin(^), where <£ is the phase difference between 
the two arms of the interferometer. 

On can write <£ as the sum of three phases: 



(p 



Bloch 



& 



kinctii 



(24) 



The Raman and kinetic phase are the same as for the usual 
Ramsey-Borde interferometer. The Bloch phase includes 
all the phases shift during the Bloch separation, due to the 
phase of the laser, the light shift of atoms in the lattice 
and the kinetic energy. 

The Bloch phase is the difference between the Bloch 
phase ^ Bloch accumulated in the path A and B. For a 
given sequence of the optical lattice (velocity and ampli- 
tude) and an initial momentum pq, we calculate the co- 
efficients c; (c.f. eq. [7]) corresponding to the evolution of 
the state \po) (resp. \po + 2hk)) using the model described 
earlier. The contribution of the first pulse to Bloch phase 
is deduced from the phases of the coefficients c~n (resp. 
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Fig. 11. Interference fringes observed using the scheme de- 
picted on Fig. [10] On this figure, the main source of noise is 
the detection noise. 



Cjv+i)- The same method is used for all Bloch pulses. This 
allow us finally to calculate the total phase for the 2 first 
beamsplitter, # { . From this phase we deduce the posi- 



tion of the wave packets (Xa/b) 



,)•!■ 



A/1 



The sensitivity 



of the interferometer can be deduced using equation f|23[) . 



4.2 Experimental results 

The experimental setup and the main results have been 
described in a previous paper |13) . We briefly remind the 
results to compare them with the theoretical model. We 
use rubidium 87 atoms, consequently the two levels cou- 
pled by the Raman transition are the \F — 2) and \F = 1) 
hyperfine level of the ground state. 

The Fig. [TTJ shows the proportion of atoms detected 
in the internal state F = 1 as a function of the frequency 
difference between selection and measurement using the 
time sequence reported on Fig. 1101 

The periodicity of the fringes (Fig. [TTJ is about 5 scp = 
50 Hz for the time sequence depicted on Fig. [TUJ 

Quantitatively, <5 sop can be compared with the period- 
icity of fringes of an interferometer without LMT Bloch 
pulses. With a delay of 5 ms between the two Ramsey 
pulses of the selection, the periodicity would have been 
200 Hz. Therefore with LMT pulses, the sensitivity of the 
interferometer is improved by about four. 

We can define an effective time T c g = l/<5 sep , which is 
of about 20 ms. Because the fringes spacing depends only 
on the separation of the two arms of the interferometer, 
the effective time T e g can be viewed as the Ramsey time 
in a regular Ramsey-Borde interferometer which leads to 
the same spatial separation of the two wave packets. 

More rigorously, T e g and the spatial separation (Xa) — 
(Xa) of the two wave packets are linked through the fol- 
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Ramsey I 



Fig. 12. Effective time T e g as a function of the total time of 
the selection. The temporal sequence described Fig. [10] is kept 
constant, except for TR amscy . Solid curve: Numerical calcula- 
tions. Inset: sensitivity gain of the LMT interferometer com- 
pared to a regular Ramsey Borde interferometer of the same 
total duration. 



lowing equation: 



r P , 



(X A ) - (X B ) ^ 
2v r 



(25) 



The dependence of T c g with the duration of the selec- 
tion sequence has been also investigated experimentally. 
The duration of the time-sequence of the interferometer 
is adjusted by changing TR amscy keeping the same ir/2- 
Bloch-Bloch sequence. The variation of T c g versus TR amsey 
is reported again on Fig. [T2] There is an excellent agree- 
ment between experiment and calculations. 

In the inset of the Fig.[12l we have also plotted the gain 
in resolution (ratio between the effective time and the real 
time) as a function of the duration of the selection. At the 
limit where TR amsC y is long, the duration of the pulses can 
be neglected and the gain will be exactly 27V + 1, i.e. 5 for 
our parameters. At shorter times, this gain is smaller. 

The experimental data fits well with the predicted time 
which is a strong evidence that the fringes observed on 
Fig. [TT] are due to the LMT pulses. 



4.3 Light shift in the interferometer 

Light shift in the interferometer are one of the main source 
of systematic effect. Furthermore they can induce a reduc- 
tion of the contrast of interference. 

The mean energy of atoms is indeed shifted when the 
lattice is switched on. This shift is not the same between 
the different bands and it depends on the intensity of the 
lattice. Our scheme (see Fig. fTUJ) is on average symmetric: 
that is the time spent by atoms in a given energy band 
is equal for each arm of the interferometer. Consequently, 
there is no associated systematic shift if the intensity seen 
by atoms is constant. However, experimentally, this con- 
dition is never realized because of temporal fluctuations 



of the laser intensity and because of motion of the atoms 
through the spatial profile of the laser beam. In the first 
case, fluctuations of the intensity will be the same for each 
atom, leading to a phase noise in the interference pattern. 
In the second case, the light shift contribution will be av- 
eraged over all the atoms, leading to a reduction of the 
fringes contrast. 

The light shift are implicitly included in the phase 
^Bioch obtained by integrating the Schrodinger equation. 
This phase depends on the intensity of the laser over time. 
For small variation of intensity, one could linearized the 
variation of the phase and compute the functional deriva- 
tive of the phase with respect to the intensity: 



o(t) = 



5<P 
W) 



(26) 



In the experiment the intensity of the laser is modulate 
by an acousto-optic modulator. The intensity is therefore 
the product of the modulation (M(t)) due to the AOM 
and the intensity without modulation Io(t). If we assume 
that the modulation for the LMT pulse does not fluctuate 
(5M(t) — 0), then the fluctuation of the phase is mainly 
due to fluctuations of Iq (t) around its constant mean value 

In the following instead of calculating -iMp- , we choose 
to calculate the functional derivative with respect to the 
normalized derivative of Io(t), 1/Iq dIo(t)/dt, that we de- 
note by <?/(£). We choose to use the derivative, because 
it is more convenient numerically and also lead to a di- 
mensionless functional derivative. At first order, one can 
therefore write that: 



5<P- 



, 9l{t) T ^ dt 



(27) 



In the following, gi is referred as the sensitivity function 
to the relative derivative of intensity. 

Numerically, gi (t) is computed by looking to the varia- 
tion of the phase <?Bioch for a small change of the intensity 
Io{t') for t' > t (Heaviside function). The temporal evo- 
lution of the function gj during the sequence described 
previously on Fig. [ID] is reported on Fig. [13] 

The inset is a zoom of gi at the beginning of the se- 
quence. In the inset, the contributions to gi of each arm 
{gf and gf) are plotted. In the experiment, atoms that are 
in excited band experience on average, the mean potential 
of the lattice, Uq/2, whereas atoms in the first band, are 
in dark region (blue detuned lattice) and less subjected to 
light shift. 

The phases in the inset of Fig. Q2] are calculated with 
cq. [5] in which the common light shift Uo(t)/2 is left out. 
This changes the values of gf and gf but not the value of 
<7j. Within this energy reference, an atom moving fast with 
respect to the lattice does not see a phase shift ((cos(2fcx)) - 
0) whereas an atom trapped into the lattice sees a phase 
shift ((cos(2fcx)) f=a —1). Therefore the main contribution 
to the phase shift originates from the lower band. Conse- 
quently, it explains why the initial increase of gi princi- 
pally arises from the first band. 
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The contrast C is given by C = e A ^ I" 1 where A<j> is the 
variance of the phase </>. Acj) is given by: 



A(j) = (Tj 



gi&t 



(29) 



where <7j, which represents the relative variance of -gl- 
over the entire atomic cloud, is calculated with 



r„\ 




(30) 



To evaluate <jj numerically, Gaussian distributions are 
used to described the position and the velocity of the 
atoms. Then <jj can be expressed as: 



a; = 4 v2er„ — f 



1 



wi 1 + ; 



(31) 



Fig. 13. Sensitivity function of the interferometer to the rela- 
tive derivative of intensity. The temporal sequence is described 
Fig. 1101 Inset: zoom on the two first LMT pulses. The contri- 
bution from the two arms of the interferometer are displayed 
separately. A constant phase shift corresponding to the aver- 
age light shift of an atom in the lattice is subtracted from both 
curves. 



This phase shift is almost compensated by the other 
arm during the second LMT pulse. Indeed, the first two 
LMT pulses do not act in a symmetric way between the 
two arms because non-accelerated atoms are not in the 
same excited band during the first and second LMT pulses. 
During the first LMT pulse, non- accelerated atoms en- 
dure a small light shift (see the variation of gf , around 
t=0.6 ms in the inset of Fig. [T21 arrow labeled a). While 
during the second LMT pulse, non-accelerated atoms are 
moving too fast with respect to the lattice to see a light 
shift (see the plot of gf which remains constant during 
the second LMT, at t around 1.2 ms, arrow labeled b). 

In the case of the interferometer with only four LMT 
pulses (see Fig. [9J), the initial phase due to the first LMT 
pulse is compensated only at the end of the interferom- 
eter, after about 10 ms, whereas in our scheme, 90% of 
this phase shift is compensate after 1 ms. The sensitiv- 
ity of our interferometer to a linear variation of the phase 
(f gi(t)dt) is then reduced by a factor of 10 compared to 
an interferometer with four LMT pulses. 

From the sensitivity function gi, the reduction of the 
contrast due to the motion of atoms during the interfer- 
ometer can be precisely calculated. For that, the phase 
shift accumulated along a trajectory must be calculated. 
At first order, it is expressed as 



1 d/ °^ 
flr at 

y /o at 



9ldt ' T ~dT 



(28) 



t=0 



in which a denotes the RMS value of the position of the 
atoms in one direction and a v the RMS value of their 
velocity, raj, represents the waist of the Gaussian laser 
beam. With ours experimental parameters (wl = 2 mm, 
a = 600 /im, a v — 2v r ), the value of Of is around 5.9 s . 

For the eight LMT pulses scheme used in our experi- 
ment, we have calculated that: J gidt sa 48 x 10~ 3 rad • s. 
Therefore, the values of A(f> and C are : Ac/) = 0.28 rad 
and C = 96%. 

With the four pulses scheme, those values become Acj> = 
2.8 rad and C — 2%. This clearly shows that in our ex- 
periment, the eight LMT pulses scheme must be used to 
observe interferences. 

We emphasize that there is a correlation between the 
velocity and the position of atoms because of the small 
time of flight before the beginning of the interferometer. 
This results in a biased intensity variation leading to a 
systematic effect in the interferometer. This effect can be 
canceled by inverting the order of the two LMT pulses for 
each beamsplitter because the sign of gi and consequently 
the one of d is reversed. 



5 Improving LMT beamsplitter 

5.1 Current limitations 

The contrast of the fringes observed Fig. [TT] is ten times 
smaller than the one without LMT pulses. This weak con- 
trast can not be explain with light shifts. In the experi- 
ment, light shift is under control when the LMT pulses are 
used on both arms of the interferometer. Other sources 
contribute to this reduction. The main one arises from 
inhomogeneities in the initial momentum distribution of 
atoms and intensity of the laser. The maximum efficiency 
computed has been obtained for an atom with a defined 
velocity and a given laser intensity. But in our set-up, the 
initial velocity distribution is spread over about 0.5 v r , 
consequently the effective efficiency for a LMT pulse is re- 
duced from 98% to 95%. Moreover, because of the initial 
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size of the cloud and the finite waist of the laser, the varia- 
tions of the intensity seen by atoms are in the order of 20%. 
As it can be easily see on Fig. [21 this effect pushes away 
the efficiency from the optimal configuration by few per- 
cents. The contribution of spontaneous emission of atoms 
trapped in the lattice is also non negligible. In our regular 
experiments with Bloch oscillations [18], atoms, located 
at the bottom of the lattice wells, are oscillating in the 
first band. As the potential of the lattice is blue detuned, 
the spontaneous emission is reduced, because the light in- 
tensity is minimal at the position of atoms |18j . Unfortu- 
nately, this demonstration does not apply to atoms in the 
excited bands. The estimated spontaneous emission rates 
of those atoms is around 4%. 

All these effects added together reduce significantly the 
efficiency of each of the 8 LMT beamsplitter used in the 
interferometer and explain the reduction of contrast that 
we have observed. However, using a laser beam with a 
larger waist, a higher intensity (and a larger detuning) 
should allow us to reduce significantly these effects. 



5.2 Sequence of a 10 recoil beamsplitter 

As mention previously, it is not possible to realize a larger 
momentum separation while keeping a good efficiency with 
the scheme of the LMT, in which constant acceleration is 
used (c.f. FigEJ. 

To overcome this limitation, a process with different 
acceleration steps can be foreseen. The first step is iden- 
tical to the one described earlier, with two Bloch oscilla- 
tions. For the second one, the depth of the lattice is raised 
to enable a larger acceleration. 

As explained previously, a long time and a relatively 
weak lattice are necessary to realize the first splitting with 
a good efficiency. For a given time, the optimal value of 
the amplitude of the lattice is a compromise between good 
adiabaticity of Bloch oscillations in the first band and a 
high Landau-Zener tunneling for the third band (n = 3). 
However, this condition is less restricting when atoms are 
in a higher band (n > 3). For those atoms, the Landau- 
Zener tunneling is higher so one can expect a better effi- 
ciency with the same duration or the same efficiency with 
a shorter duration. That is, a larger acceleration can be 
use to increase separation with a good efficiency. 

The realization of a beam splitter with ten oscillations 
is presented on Fig. 1141 The initial separation with two 
BOs is achieved in 130 fis while the second separation 
with eight BO is done in only 160 /us. The maximal ef- 
ficiency is calculated to be about 98%. This is very close 
to the efficiency with two BOs beamsplitter. The main 
losses arise from the initial splitting of two recoils. The 
additional losses come from the lack of adiabaticity dur- 
ing the change of the lattice intensity. As Bloch oscillations 
is a very efficient process, no significant losses arises from 
the additional recoils transferred after the first two oscilla- 
tions. The main limitations will come from the reduction 
of contrast due to intensity fluctuations. 

To estimate this limitation, the sensitivity function gz, 
with ten BOs, is plotted on Fig.[T5] It is significantly larger 
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Fig. 14. Scheme of a beamsplitter with 10 oscillations. 
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Fig. 15. Sensitivity function gi of the interferometer using a 
10 Bloch oscillations LMT pulse, as described in Fig. 1141 



than the one with two BOs because of the increase of the 
depth of the lattice. The sensitivity to a linear variation of 
intensity J gjdt is about 95 x 10~ 3 rad • s. That is twice as 
much as the one for the interferometer with two BO (see 
section I3TTJ) . This is the reason why the improved LMT 
was not implemented yet on our current apparatus (the 
visibility of fringes for the "two BOs scheme" is already 
quite weak). 



6 Conclusion 

In this paper, we have theoretically and experimentally 
studied the implementation of a large momentum beam- 
splitter in a Ramsey-Borde interferometer. We have real- 
ized a separation of 10 recoils between the two arms of the 
interferometer. By reducing the contrast of the fringes, 
light shift is the main limitation of this interferometer. 
Nevertheless, we have developed a method to reduce this 
effect by a factor of ten and so realize an interferome- 
ter with a separation of ten recoils between the two arms. 
Our model is in very good agreement with the experiment, 
in which a gain of four in the resolution, compared to a 
usual interferometer, was observed. This method seems to 
be very promising for the realization of an interferometer 
with a separation of several tens of recoil velocities. 

The limitation due to light shift has been strongly re- 
duced by accelerating successively both arms of the in- 
terferometer. Moreover a way to systematically cancel it 
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would be to accelerate both arms of the interferometer si- 14. 
rrmltaneously (instead of doing it successively). This can 
be done by applying two counter-propagating accelerated 15. 
lattices. This scheme is already used by the Stanford/Berkeley^, 
group [5]. In their experiment, a high order Bragg beam- 
splitter is initially used, before the double BO accelera- 17 - 
tion. We are also currently investigating this possibility of 
using the double BO acceleration but with an initial ve- 
locity separation given by a two photon Raman pulse. It 
seems that only two recoils is two small to perform a good 
separation and a higher initial separation may be required. 
A relevant method could be to use the double-diffraction 
technique described in Ref. [10] which allow an initial sep- 
aration of 4 recoils. In this scheme, light shift would be 
completely suppressed. One of the remaining source of 
noise will be then the phase/vibration noise. This noise 
can be calculated with a numerical model identical to the 
one used in this paper. 
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